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Abstract 

A detailed compilation of uncertainties in the MS bottom quark mass m^m;,) 
obtained from low-n spectral sum rules at order is given including charm 
mass effects and secondary bb production. The experimental continuum re- 
gion above 11.1 GeV is treated conservatively. An inconsistency of the PDG 
averages for the electronic partial widths of T(4S) and T(5S) is pointed out. 
From our analysis we obtain m^(jn^) = 4.20 ±0.09 GeV. The impact of future 
CLEO data is discussed. 



Introduction 



Present data from B factories on inclusive decays already require precise knowledge 
of the bottom quark mass parameter with a numerical precision of the order 50 MeV 
with a reliable estimate of the uncertainty. This will become even more acute in the 
future when more data becomes available. Using methods based on perturbative 
QCD, there have been several approaches in the past aiming at uncertainties of 
less than 100 MeV. The most frequently used method is based on large-n (n > 4) 
moments of the bb production cross section in e + e~ annihilation 

Pn- j^RM, (1) 

where R b i = a(e + e~ — > bb + X)/cr(e + e~ — > u + u~) and the contributions from the 
virtual Z are neglected. From fitting moments obtained from experimental data 
to the corresponding theoretical expressions the bottom mass can be determined in 
threshold schemes (for recent reviews see Refs. Large- n ( "non-relativistic" ) mo- 
ments have the advantage that the badly known experimental continuum bb cross 
section above the T(6S) is strongly suppressed in comparison to the rather well 
known resonance region. However, the order a 2 s corrections in the framework of the 
non-relativistic expansion are generally quite large and small uncertainties below 
100 MeV can only be achieved with additional assumptions on higher order correc- 
tions. This indicates that an improvement in the treatment of the non-relativistic 
bottom quark dynamics might be needed. 

A different method uses low-n ( "relativistic" ) moments [§, |6| where n < 4. In 
the recent past, relativistic moments have been used less frequently because the 
badly known continuum region represents a major source of uncertainty that is not 
reducible without additional assumptions. Theoretically, the usual loop expansion 
in powers of a s can be employed since for small n the bottom dynamics is relativis- 
tic. Here, the MS mass is an appropriate mass definition to be used and extracted. 
In contrast to large-n moments, the low-n moments show a quite good perturbative 
behavior. A recent analysis by Kiihn and Steinhauser used moments at order 
a 2 s . Adopting a theory-driven perspective, the experimental continuum data above 
the T(6S) was obtained from theoretical results for R bb , basically eliminating uncer- 
tainties from the continuum region. In Ref. [§] the important conclusion was drawn 
that, using the strategy of Ref. 0, a substantially more accurate measurement of 
the T(1S) - T(6S) region at CLEO |§ could result in an uncertainty in m b (m b ) of 
only 30 MeV. 

It is the main purpose of this paper to give a detailed compilation of all sources 
of uncertainties in the bottom MS mass m b (m b ) obtained from low-n moments, 
including a more conservative treatment of the experimental continuum region. We 

1 Recently, the bottom mass has been determined in different threshold schemes from mo- 
ments of inclusive semileptonic and radiative B meson decay spectra with an uncertainty of about 
100 MeV, which is dominated by experimental errors. || (See also Ref. |Q.) 
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believe that this compilation can contribute to a more differentiated view on the 
current uncertainties in m^mj) from low-n moments and on the impact of new 
more precise data in the T resonance region from CLEO. In our analysis we also 
include the contributions from secondary bb production from gluon splitting and the 
effects of the non-zero charm quark mass, which have to our knowledge not been 
taken into account before. Both effects turn out to be small. Finally, we point out an 
inconsistency in the way the PDG has treated the original results for the electronic 
partial widths of T(4S) and T(5S) from CUSB and CLEO JTTJ, which leads to a 



contribution to the experimental moments P n that is smaller than the contributions 
obtained from the original data in that energy region, both from CUSB and CLEO. 



Theoretical Moments 

For the QCD parameters used in this work we adopt the MS renormalization scheme 
and the convention that the bottom quark participates in the running (n/ = 5). The 
masses of the quarks in the first two generations are set to zero. In terms of rfib(^) 
the moments in the OPE, including the known perturbative corrections to order a 2 
and the contribution from the dimension four gluon condensate, take the form 
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where r = m c /mb. In terms of the more specific choice of m^rrib), the moments 
have the simpler form 
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where /?o = 11 — 2/3rif. The Born and order a s terms of the moments are known 
since a long time 0, |6|] and the order a 2 s contributions for primary bb production 



for massless light quarks have been determined in Ref. ||12| . We have cross-checked 
these contributions with the explicit expressions for the corresponding contributions 
to R bb given in Ref. [13| . The order a 2 s contributions to the moments from secondary 
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Table 1: Coefficients of the theoretical expressions for the moments P„ to order a 2 
for massless light quarks. 
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-0.0466 
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-0.0109 
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Table 2: Corrections due to the non-zero charm quark mass to the order a 2 coefficient 
f 20 for r = me/nib. 



bb production, where the bb pair is produced through gluon radiation off light quarks, 



has been computed from the corresponding results for the i?-ratio given in Refs. |L4] , 
15|. These contributions only affect the coefficient f 20 . The coefficients of the gluon 



condensate have been taken from Ref. [|l^, [T^J . For convenience, the numerical results 
for the coefficients in Eqs. (0) and (|3]) for n = 1,2,3,4 and m c = (r = 0) are 
collected in Tab.|l| The numbers for the / n 's agree with Ref. [0] up to a convention 
dependent factor of 1/4, except for the results for /^°(0), which are slightly larger 
accounting for the contributions from secondary bb production. The effects of the 
non-zero charm quark mass are generated either through virtual gluon self-energy 
effects or through real primary or secondary associated charm production. The 
corresponding contributions to R h i for arbitrary mass constellations have been given 
in Refs. [|i~3| , |i~4"f . Numerical values of the charm quark mass corrections to the 
coefficient f%° are displayed in Tab.|| for values of r between 0.1 and 0.5. We note 
that the numbers given in Tab.|2| also include the non-zero charm mass effects in the 
bottom quark pole-MS mass relation TJ| . 
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T(1S) 


0.766(29) 


0.856(32) 
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2.467(26) 
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Table 3: Individual contributions to the experimental moments including uncertain- 
ties. The contribution from a resonance k has been determined in the narrow width 
approximation, (P n )k = 97r T| +e /[a (10 GeV) M^ n+1 ], where for the electromagnetic 
coupling [a (10 GeV)] -1 = 131.8 has been adopted. 



Experimental Moments 

For the contributions to the experimental moments from the T(1S), T(2S), T(3S) 
and T(6S) we use the averages for masses and e + e~ widths given by the PDG flTBf . 
In Tab. |] a collection of all contributions to the moments including uncertainties is 
given. The averages for the T(1S), T(2S) and T(3S) are dominated by data from 



ARGUS (19| and the averages for the T(6S) are from results from CUSB ||T0[ and 



CLEO [11 



In the 4S-5S region between 10.5 and 10.95 GeV there have been measurements 



from CUSB [10| and CLEO [|Tl|]. We find it remarkable that both experiments 
observed an additional resonance-like enhancement between the T(4S) and T(5S) 
at about 10.7 GeV. Whereas CUSB fitted for T(4S) and T(5S) resonances within in 
a coupled channel model, the CLEO experiment was fitting an additional resonance, 
called u B* n , at m B * = 10.684 ± 0.013 GeV with an e+e" width of Y e ^ e ~ = 0.20 ± 
0.11 keV. As a consequence the e + e~ widths for T(4S) and T(5S) obtained from 
CLEO were systematically smaller than those from CUSB. In the PDG compilation, 
however, the existence of the U B*" contribution in the CLEO analysis was ignored 
when the averages for T(4S) and T(5S) have been determined.^ As a consequence, 
the contribution to the moments from the region between 10.5 and 10.95 GeV using 



2 The PDG number for the electronic width of the T(4S) is a weighted average of the results 
from CUSB Jh}, CLEO || and ARGUS @, where the ARGUS result agrees better with the one 
from CUSB. 
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the PDG averages is, although compatible within errors, systematically lower than 
the contributions one obtains using the numbers given in the original CLEO and 
CUSB publications (see Tab.|3|). For our analysis we decided to ignore the PDG 
averages and to take the average of the original CLEO (including the ll B* n ) and 
the CUSB contributions to the moments for this energy region. We adopted the 
larger CLEO error assuming that the respective uncertainties are correlated. We 
believe that this conservative treatment is justified as long as the situation is not 
clarified. We note that the PDG treatment of the 4S-5S region does not affect the 
results for the hadronic vacuum polarization effects of aiMz) and (g — 2)^ because 
the corresponding differences are much smaller than the total uncertainties. For the 
T(6S) we used the PDG averages since it is unlikely that the different treatment of 
the enhancement observed between T(4S) and T(5S) has affected the fits above the 
T(5S). 

There is no direct experimental data for a(e + e~ — > 7* —>■ bb + X) in the region 
above 11.1 GeV. However, there are measurements of the total hadronic cross section 
taken by a number of experiments up to energies close to Mz that are compatible 
with the Standard Model predictions. Furthermore, from measurements of R b at 
the Z pole by LEP and in the region between 133 and 207 GeV by LEP20, it is 
known that perturbative QCD agrees with the data for bb production to about 1% 



at Mz and to about 10% in the LEP2 region [20fl. It is therefore not unreasonable 



to estimate the experimental contribution to the moments from above the T(6S) 
from perturbation theory itself. Although order a\ corrections to R b i in the high 
energy expansion are known, we use the perturbative bb cross section to order a 2 to 
estimate the continuum contributions because the theoretical moments are likewise 
only available to order a 2 . A much more subtle question is how to estimate the "ex- 
perimental" uncertainties in this region. The approach of Ref. assumes that the 
experimental data for the bb cross section including errors lie within the theoretical 
predictions and uses the small theoretical errors. This approach is quite similar to 
using finite energy sum rules |21j where an upper "duality" cutoff s max > 11.1 GeV 
is used in the integral in Eq. (fl|). 

In Tab. [3] we have displayed the continuum contributions to the moments to 
order a 2 . The charm quark mass has been set to zero. We have subdivided the 
continuum contribution into three parts coming from 11.1 — 12.0 GeV (region 1), 
12 GeV — Mz (region 2) and Mz — 00 (region 3) in order to visualize the impact of 
the various energy regions. The theoretical errors shown in Tab. |3] come from varying 
the strong coupling in the range a s (Mz) = 0.118 ± 0.003 and the MS bottom mass 
in the conservative PDG range m&(ra&) = 4.2 ± 0.2. The renormalization scale 
\i has been varied between 2.5 and 10 GeV. For the running of a s and m& four- 
loop renormalization group equations have been used. Our theoretical errors are 
larger than in Ref. where also order contributions have been included [22 



3 At the Z pole Rb is defined as the ratio of the total b quark partial width of the Z to its total 
hadronic partial width, T^/Thad, and for LEP2 energies it is defined as the ratio of the total bb 
cross section to the total hadronic one. 
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Our central values have been obtained from the average of the respective extremal 
values. Region 1 has been displayed separately because data for the bb cross section 
could potentially be collected there by CLEO || in the near future. 

Since the continuum region above 11.1 GeV is unsuppressed and constitutes a 
sizeable contribution to the experimental moments^, using the small theory errors 
leads to a considerable model-dependence of the bottom quark mass. For our final 
error estimate of m b (m b ) we adopt a 10% correlated (relative) error for the continuum 
regions 2 and 3 and ignore the respective theory errors given in Tab.|3|. This choice 
is, in principle, as arbitrary as using the theoretical errors (or no errors at all), but 
should reduce the model-dependence to an acceptable level. For region 1 we use 
the theory error since here the variation of the MS mass in the conservative PDG 
bounds has the largest impact and leads to a variation of more than 10%. 

Uncertainties in ra&(ra&) 

For the determination of m b (rrib) and the uncertainties we have used 4 methods: 

1. The bottom mass m b (m b ) is determined from single moment fits (n = 1, 2, 3) 
using Eq. (|]). 

2. The bottom mass m b (fi) is determined from single moment fits (n = 1, 2, 3) us- 
ing Eq. (0) and m b (m b ) is computed subsequently using renormalization group 
equations. 

3. The bottom mass m b (m b ) is determined from fits to ratios P n /P n+ \ (n = 1, 2) 
using Eq. (|]). 

4. The bottom mass m b (/i) is determined from fits to ratios P n j P n +i (n = 1, 2) us- 
ing Eq. (g) and m b (m b ) is computed subsequently using renormalization group 
equations. 

For the analysis we employed only moments for n = 1, 2, 3 to avoid the large higher 
order contributions ~ (a s ^/n) k that are characteristic for the large-n moments and 
need to be summed. For method 3 and 4 we did not expand the perturbative series 
in the theoretical ratios P n /P n+ i. We checked that expanding the theoretical ratios 
has only very small effects on the results. We employed four-loop renormalization 
group equations and used a s (M z ) = 0.118 ± 0.003, m c = 1.3 ± 0.2 GeV, (^G 2 ) = 
(0.024 ± 0.024) GeV 4 as theoretical input. The renormalization scale /i was varied 
between 2.5 and 10 GeV. 

In Tab. | the results of our analysis for m b (Jn b ) are displayed in detail. The 
table shows the respective central values (in units MeV), which were obtained using 

4 We note that the relative contribution to the low-n cc spectral moments coming from contin- 
uum energies above 4.5 GeV is considerably smaller than for the continuum region above 11.1 GeV 
in bb spectral sum rules, see e. g. Ref. M. 
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Table 4: Central values and uncertainties for m h {m b ) in units of MeV based on the 
methods described in the text. 



the experimental and theoretical central values given before and /i = 5 GeV. The 
central values obtained with the four methods are within 15 MeV around 4.20 GeV. 
All uncertainties (in units of MeV) are presented separately. The experimental 
errors correspond to the uncertainties given in Tab.|3], where for the 4S-5S region 
our conservative CUSB-CLEO average has been used. The uncertainty from the 
4S-5S region constitutes the largest experimental error from the resonance region. 
The errors from the continuum regions indicated with a subscript "th" are obtained 
from the corresponding theory errors shown Tab.|3|. For the continuum regions 2 
and 3 also the errors coming from a 10% deviation from the theory prediction are 
displayed (having no subscript). We note that the latter errors scale roughly linearly, 
i.e. assuming a 5% (20%) deviation the error decreases (increases) by a factor of two, 
etc.. This illustrates how strongly the bottom quark mass depends on assumptions 
for the experimentally unknown bb continuum cross section above the T(6S). The 
theoretical errors have been obtained by varying each of the theoretical parameters 
in the ranges given above while the respective other parameters were fixed to their 
central values. For method 3 and 4, which are based on fitting ratios P n / P n+ i, the 
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same theoretical input parameters have been used for the moments in the numerator 
and those in the denominator. If a s (Mz) is chosen independently for the moments 
in the numerator and denominator, the errors for method 3 (4) are 38 (26) MeV for 
n = 1 and 39 (13) MeV for n — 2. If /x is chosen independently for the moments 
in the numerator and denominator, the errors for method 3 are 52 MeV for n = 1 
and 63 MeV for n = 2. We found that, except for the variations of fi, all resulting 
errors scale linearly with changes of the input parameters. The errors coming from 
variations of fi have been chosen to be the larger ones of the two deviations obtained 
in the ranges 2.5 GeV < /i < 5 GeV and 5 GeV < fi < 10 GeV. Note that the overall 
shift in the central value of m^m;,) coming from the gluon condensate contribution 
is between —0.1 and —1 MeV. The shift caused by the non-zero charm quark mass 
is between —1 and —3 MeV, which is an order of magnitude smaller than for large-n 
moments, where the charm mass effects are enhanced by a factor l/a 2 s ~ n [23 



Using the PDG average for the 4S-5S region instead of the CUSB-CLEO average 
(see Tab.|3|), m b {m b ) is shifted by 10 to 15 MeV for method 1 and 2 and by 2 to 
6 MeV for method 3 and 4. We consider the numbers given in Tab.|] as the main 
result of this work. 

In order to obtain combined errors from the uncertainties of the resonance data 
we treated one half of each error as correlated (being added linearly) and the other 
half uncorrelated (being added quadratically) because all data came from e + e~ ma- 
chines with common systematic uncertainties and, roughly, systematic and statisti- 
cal uncertainties were found to have comparable sizes [fTl] . In Tab.^ the resulting 
combined resonance errors are shown in the line below the numbers for the T(6S). 
The theoretical uncertainties including the errors adopted for the three continuum 
regions have been combined linearly since they do not have any statistical meaning 
and, in particular, the division into three continuum regions is completely arbitrary. 
Note that for regions 2 and 3 the uncertainties from the 10% variation of the theory 
prediction have been adopted and not the smaller theoretical errors. The resulting 
combined error is displayed in the line below the numbers for Sfi. To obtain our total 
error (last line in Tab.fD we added the resonance, the continuum and the theory 
errors linearly. As expected we find that the uncertainties from the continuum have 
the largest impact for the moments with n = 1 and 2 and that they are partially 
canceled when ratios are used for the fitting. However, the theoretical errors are 
larger for fits with ratios than with single moments, if the theoretical parameters 
are chosen independently for numerator and denominator. In general, the total error 
decreases for larger n. This trend does, however, not continue for higher values n > 3 
particularly for methods 2 and 4. Compared to the results of Ref. J7| our total errors 
are much larger, particularly for fits involving Pi and Pi- This is mainly because 
we adopted more conservative errors for the continuum region in the experimental 
moments, and we combined uncertainties linearly, when they cannot be treated sta- 
tistically. Since we believe that P3 can be computed reliably using Eqs. (0) and (^), 
we adopt the average of the total errors in the third column of Tab. as our final 



8 



estimate for the uncertainty and obtain (rounded to units of 10 MeV) 

m b (m b ) = 4.20 ± 0.09 GeV. (4) 

We do not take into account the small 50 MeV error from P2/P3 for method 3 
because the small error from variations of /i only persists if the same /1 is chosen for 
both moments. For an independent choice the error is considerably larger (see text 
above). If a 20% relative uncertainty is assumed for the continuum regions 1, 2 and 
3, the final error increases by 20 MeV. Our result in Eq. (f|) is compatible with the 
result from Ref. 0. Our result is also compatible with results from fitting large- 
n moments at NNLO in the non-relativistic expansion and with recent results, 
using different methods, obtained in Refs. |2T|, . 

Based on the numbers given in Tab.^| and f| it is straightforward to discuss the 
impact of improved measurements of the resonance parameters at CLEO. Assuming 
improved measurements for the electronic widths of T(1S), T(2S) and T(3S) |§ at 
the level of 2%, the combined resonance errors shown in Tab. ||| would be reduced by 
a factor 2/3, which would reduce the total error (obtained from the third column) by 
about 10 MeV. An improved measurement of the 4S-5S region and the T(6S) with 
the same precision would result approximately in a reduction of the total error by an 
additional 10 MeV. A further reduction of the error below about 70 MeV, however, 
will be difficult to achieve without real experimental data for a(e + e~ —>■ bb) in the 
continuum region above the T(6S) with a precision of better than ten percent. 

Conclusion 

We have given a detailed analysis of the uncertainties in the MS bottom quark 
mass m b {m b ) using low-n spectral moments of the cross section cr(e + e~ — > 7* — > 
bb). For the experimental moments we employed experimental data for the T(1S) - 
T(6S) resonances and the order a 2 QCD predictions for the continuum region above 
11.1 GeV. For the 4S-5S region between 10.5 and 10.95 GeV we found that the 
PDG averages for the electronic partial widths of T(4S) and T(5S) based on data 
from CUSB, CLEO and ARGUS contain an inconsistency, stemming from the fact 
CLEO also assumed the existence of an additional resonance at about 10.7 GeV 
which was ignored in the averaging procedure. For our analysis we therefore used 
the original CUSB and CLEO results. For the continuum region above the T(6S) 
we assumed a 10% error in our final error estimate. For the theoretical moments we 
used perturbative results at order a 2 including also the contributions from secondary 
bb production and finite charm mass effects. As our final result we get rn^im^) = 
4.20 ± 0.09 GeV and we conclude that more precise data for the electronic partial 
widths of the T resonances at CLEO could reduce the error by about 20 MeV. 
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